Belle Preprint 2004-15 
KEK Preprint 2004-19 

Measurement of Branching Fraction and CP Asymmetry in B^ -^ p+7r° 
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We report a measurement of the branching fraction for the decay B"*" -^ p^-K^ based on a 140 fb~^ 
data sample collected with the Belle detector at the KEKB asymmetric e^e" collider. We measure 
the branching fraction B(B+ -^ /o+tt") = (13.2 ± 2.3(stat.)li j(sys.)) x 10"^, and the CP-violating 



asymmetry Acp{B^ 



T^ON 



p^-K 



= 0.06±0.19(stat.) 



lH^(sys.). 
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Recent precise measurements of sin 2(j)i [l|, I3] confirm 
the prediction of the Kobayashi-Maskawa model |3] for 
CP violation. It is of great importance to test this theory 
further with complementary measurements, such as those 
of the other unitarity triangle angles 02 and cf)^ |j| . 

At the quark level, the decays B -^ pir occur via & — > u 
tree diagrams and can be used to measure (/>2. How- 
ever, because of the presence oi b ^ d penguin (loop) 
diagrams, a model independent extraction of 02 from 
time-dependent CP-asymmetry measurements requires 
an isospin analysis of the decay rates of all the gn de- 
cay modes 5i]. The decay channels B^ -^ p^ir^ fg and 
_B° -^ p^TT^ have already been measured J7|. Evidence 
for the _B" -^ p^-K^ mode, which is expected to be small, 
has been reported by Belle [g] with a rate higher than an 
upper bound from Babar [9j . The remaining decay mode, 
5+ -^ P^tt", has two neutral pions in the final state that 
make its measurement an experimental challenge. Re- 
cently, the BaBar group reported the observation of this 
mode 0. 

In this Letter, we report measurements of the branch- 
ing fraction and the CP-violating charge asymmetry for 
the P+ -^ p+7r° decay mode. The results are based on 
a 140 fb^^ data sample containing 152.0 x 10^ B me- 
son pairs collected with the Belle detector at the KEKB 
asymmetric-energy e+e^ collider [lOJ operating at the 
T(45') resonance {^/s = 10.58 GeV). The production 
rates of B'^B~ and B^B^ pairs are assumed to be equal. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a three-layer silicon vertex de- 
tector, a 50-layer central drift chamber (CDC), an array 
of aerogel threshold Cerenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of 
CsI(Tl) crystals (ECL) located inside a superconducting 
solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to 



detect Kl mesons and to identify muons . The detector 
is described in detail elsewhere [ll|. 

For charged pion and kaon identification, specific ion- 
ization measurements (dE/dx) from the CDC are com- 
bined with the responses of the ACC and TOF systems 
to form likelihoods Ljr and Lk- We distinguish pions 
from kaons by applying selection requirements on the 
likelihood ratio, L^/(L7r + Lk)- Similarly, electrons are 
identified by means of a likelihood based on ECL mea- 
surements, dE/dx information from the CDC, and the 
responses of the ACC. 

The final state for the signal consists of a charged pion 
track and two tt" -^ 77 candidates. The charged track 
is required to have a transverse momentum px > 0.1 
GeV /c and to be consistent with an origin within 0.1 cm 
in the radial direction and 5 cm along the beam direc- 
tion of the interaction point (IP). In addition, the charged 
track is required to be positively identified as a pion, and 
not be consistent with the electron hypothesis. Candi- 
date 7r° mesons are reconstructed from pairs of photons 
that have an invariant mass within ±3cr of the nomi- 
nal tt" mass, where the photons are assumed to origi- 
nate from the IP, and the tt*^ resolution a varies in the 
range 5.3 MeV - 7.0 MeV depending on its momentum. 
The energy of each photon in the laboratory frame is 
required to be greater than 50 MeV for the ECL bar- 
rel region (32° < < 129°) and 100 MeV for the ECL 
endcap regions (17° < 61 < 32° or 129° < 6* < 150°), 
where 9 denotes the polar angle of the photon with re- 
spect to the beam line. The 7r° candidates are kine- 
matically constrained to the nominal tt'^ mass. In or- 
der to reduce the combinatorial background, we only ac- 
cept 7r° candidates with momenta p^^o > 0.35 GeV/c 
in the e'*"e~ center-of-mass system (CMS). We select 
p+ — > TT+TT^ decay candidates with invariant masses in 
the range 0.62 GeV/c^ < M(7r+7rO) < 0.92 GeV/c^. 

P+ -^ P+tt" candidates are identified using the beam- 



constrained mass Mbc = V-E-bcam ~ Pb^ ^^'^ ^^^ energy 
difference AE = Eb — ii^beam, where -Ebeam is the CMS 
beam energy, and pB and Eb are the CMS momen- 
tum and energy, respectively, of the B candidate. The 
AE distribution has a tail on the lower side caused by 
incomplete longitudinal containment of electromagnetic 
showers in the Csl crystals. We accept events in the re- 
gion Mbc > 5.2 GeV/c2, -0.4 GeV < A£; < 0.4 GeV 
and define a signal region as 5.27 GcV/c^ < A/bc < ulate the negative AE region, and B^ 



select the candidate with the largest Fisher discriminant. 
The MC-determined efficiency with all selection criteria 
imposed is found to be 4.4%. 

Backgrounds from B decays are investigated with MC 
simulation. For & — > c decay processes, no signal-like 
peak is found in either the Mbc or AE distribution. 
Among the much rarer charmless decays, the dominant 
backgrounds are from B^ — > p~^ p~ decays, which pop- 



5.29 GeV/c^ and -0.20 GeV < AE < 0.07 GeV. 

The continuum process e+e^ ^ 99 (<? = u,d,s,c) is 
the main source of background and must be strongly sup- 
pressed. One method of discriminating the signal from 
the background is based on the event topology, which 
tends to be isotropic for BB events and jet-like for qq 
events. Another discriminating characteristic is 0b-, the 
CMS polar angle of the B flight direction. B mesons are 
produced with a 1 — cos^ 9b distribution while continuum 
background events tend to be uniform in cos 6b ■ We re- 
quire I cos^thrl < 0.8, where ^thr is the angle between the 
thrust axis of the candidate tracks plus neutrals and that 
of the remaining tracks in the event. This distribution 
is flat for signal events and peaked at cos6'thr = ±1 for 
continuum backgrounds. We use Monte Carlo (MC) sim- 
ulated signal and continuum events to form a Fisher dis- 
criminant based on modified Fox- Wolfram moments [l2| 
that are verified to be uncorrelated with Mbc, AE and 
variables considered later in the analysis. Probability 
density functions (PDFs) derived from the Fisher dis- 
criminant and the cos Ob distributions are multiplied to 
form likelihood functions for signal (Cg) and continuum 
(Cqq); these are combined into a likelihood ratio TZs = 
Cs/{Cs + J~-qq)- Additional discrimination is provided by 
the 6-flavor tagging parameter r, which ranges from to 
1 and is a measure of the likelihood that the b flavor of the 
accompanying B meson is correctly assigned by the Belle 
flavor-tagging algorithm [ij. Events with high values of 
r are well-tagged and are less likely to originate from 
continuum production. We define a multi-dimensional 
likelihood ratio MDLR = Cf^^^/{Cf°^^ + 'C*!^^^), 
where Cf^^^^ denotes the likelihood determined by the 
r-TZs distribution for signal and £^^^^ is that for the 
continuum background. We achieve continuum suppres- 
sion by requiring MDLR > 0.9. In this way we reject 
99% of the continuum background while retaining 30% 
of the _B+ -^ p+tt" signal. 

Since B -^ p+Tr" is a pseudoscalar -^ vector + 
pseudoscalar process, the p hehcity angle 0hci, defined 
as the angle between an axis anti-parallel to the B fiight 
direction and the tt^ flight direction in the p rest frame, 
has a cos^^hei distribution. We require |cos6'hoi| > 0.3 
for further background suppression. 

In the M\^c-AE signal region, about 8% of the events 
have multiple candidates. We choose the candidate that 
has the minimum sum of x^ for the mass constrained 
TT*^ fits; in cases where candidates have the same x^, we 



7r"7r" de- 
cays, which populate the positive AE region. Monte 
Carlo studies indicate that potential backgrounds from 
_B+ -^ o-fn'^, af -^ p+tt" have AE and Mbc distributions 
similar to those for S" -^ P~^P~ decays and are accounted 
for by the latter component of the fit. In addition, the 
contamination from other possible rare B decays is taken 
into account in the systematic error. 

We extract the B+ -^ p+7r° signal yield by apply- 
ing an extended unbinncd maximum-likelihood fit to the 
two-dimensional Mbc-AiJ distribution. The fit includes 
components for signal plus backgrounds from continuum 
events, 6 — > c decays, -B" -^ P^P^ and i?° -^ tt^tt". 
The PDFs for signal, B° -^ p+p' and B° -^ 7r°7r° are 
modeled by smoothed two-dimensional histograms ob- 
tained from large MC samples. The signal PDF is ad- 
justed to account for small differences observed between 
data and MC for high-statistics modes containing tt^'s, 
i.e., B+ -^ W{K+tt-tt°)tt+ for A/bc, and D" -^ T:^n° 
for AE, where we use tt*^ mesons in similar momentum 
ranges to those from B'^ — > p+Tr*^ decay. The continuum 
PDF is described by a product of a threshold (ARGUS) 
function 13] for Mbc and a first-order polynomial for 
AE, with shape parameters allowed to be free. Back- 
ground from generic 6 — > c decays is represented by an 
ARGUS function for Mbc and a third-order polynomial 
for AE with shape parameters determined from MC. In 
the fit, all normalizations are allowed to float, except for 
the tt^tt" component, which is flxed at a MC-determined 
value based on recent Belle [13 and BaBar |lj| measure- 
ments. 

Figure n shows the final event sample and fit results. 
The six-parameter (four normalizations plus two shape 
parameters for continuum) fit gives a signal yield of 
87 ± 15 events. The statistical significance of the signal, 
defined as y'— 2 ln(>Co/'Cmax), where £max is the likeli- 
hood value at the best-fit signal yield and Cq is the value 
with the signal yield set to zero, is 8.1cr. The level of 
the B^ -^ P^P~ background determined from the fit is 
in good agreement with MC-expectations based on a re- 
cent measurement of the branching fraction for this decay 
mode K3\. 

To verify that the signal we observe is due to 5+ -^ 
p+TT^ decay, we examine the helicity and M(7r+7r°) dis- 
tributions. Figure El shows the hehcity angle distribution 
for signal yields determined from Mbc-A£^ fits, which is 
consistent with that for signal MC events. The distribu- 
tion for continuum is approximately flat. 
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FIG. 1: The upper plot is the Mbc projection for events in 
the AE signal region -0.2 GeV < AE < 0.07 GeV; the lower 
plot is the AE projection for events in the Mbc signal region 
5.27 GeV/c^ < A/bc < 5.29 GeV/c^ The solid curve shows 
the results of the fit. The signal component is shown as a 
dashed line. The continuum background is shown as a dotted 
line. The sum of 6 — > c and continuum components is shown 
as a dot-dashed line. The hatched (dark) histogram represents 
the B° -^ p+p- {B° -^ Tv°n°) background. 
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FIG. 2: Data points show the background-subtracted cos^hci 
distribution for the p; hatched histogram shows the distribu- 
tion for signal MC. The asymmetry in these distributions is 
due to the 7r° momentum requirement. 



Figure 121 shows the signal yields extracted from Mbc- 
AE fits appHed to individual M(7r+7r°) bins. A x^ 
fit to the background subtracted M(7r+7r'') distribution 
is performed with a p plus a non-resonant tttt com- 
ponent included. The p component is represented by 
a Breit-Wigner function with mass and width fixed at 
their known values [la|. The non-resonant tttt com- 
ponent is described by a second-order polynomial with 
shape parameters determined from i?+ -^ tt+tt^tt*' MC 



events, where the final state particles are distributed uni- 
formly over phase space. The fit gives the fraction of 
non-resonant decays in the 0.62 GeV/c^ < M(7r+7r") < 
0.92 GeV/c^ p signal region as (5.8 ± 4.8)%. The non- 
resonant yield increased by la is treated as a systematic 
uncertainty. 
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FIG. 3: Data points show the background-subtracted 
A4{Tr'^n'^) distribution, the dashed (dot-dashed) line is for p 
signal (non-resonant) component of the fit, the solid line is 
their sum. 



We consider systematic errors in the branching fraction 
of the decay 5+ -^ p+7r° that are caused by uncertain- 
ties in the efficiencies of track finding, particle identifi- 
cation, 7r° reconstruction, continuum suppression, fitting 
and the possible contribution from non- resonant decays. 
We assign a 1.2% error for the uncertainty in the track- 
ing efficiency. This uncertainty is obtained from a study 
of partially reconstructed D* decays. We also assign a 
0.8% error for the particle identification efficiency that is 
based on a study of kinematically selected £)*+ -^ D'^tt^, 
D^ -^ K~TT~^ decays; an 8.0% systematic error for the 
uncertainty in the two-7r° detection efficiency that is de- 
termined from data-MC comparisons of 77 ^ tt'^tt^tt^ with 
r] —> TT+TT^Tr" and 77 — > 77; a 5.1% systematic error for 
continuum suppression that is estimated from a study 



of S+ 



and D'^ -^ K^n 7r° decays; a system- 



atic error of Iq 4^ that is obtained from changes in sig- 
nal yields that occur when each parameter of the fitting 



functions is varied by ±1(t; 
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-10.5 



% systematic error 



to account for a possible contribution from non-resonant 
decays. Moreover, a 1% error due to backgrounds from 
charmless B decays other than B ^ pp and i? — > tt'^tt'^ is 
estimated by fitting the data with an additional compo- 
nent with the yield fixed at the MC-expected value. The 
change in signal yield is taken as a systematic error to 
account for this contamination. We also include a 0.5% 
error for the uncertainty in the number of BB events in 
the data sample. We obtain the branching fraction 



6(B^ 



p+TT°) = (13.2±2.3(stat.)ti-j(sys.)) x 10" 



Direct CP violation would be indicated by an asymme- 
try in the partial rates for B~ -^ p~n'^ and B"*" 



p+^°: 



A, 



CP 



T{B- 



p-7r") - ^(B^ 



-tt") 



T{B- 



p-7rO) + r(B+ ^p+ttO)' 



The B^ -^ p^TT^ candidates are self-tagged, but this tag- 
ging is provided by a single low-momentum charged pion 
that has a relatively large combinatorial background. 
This produces some wrong-tagging, which results in a 
dilution of the true asymmetry: ^obs = (1 ^ '^w)Acp, 
where .4obs is the observed asymmetry and w is the frac- 
tion of events that are incorrectly tagged, determined 
from MC to be w = 0.091 ± 0.015. 

We perform a simultaneous fit to extract the charge 
asymmetry by using the same components as shown in 
Fig. ^ and introducing asymmetry parameters into the 
fit for the _B — > pvr signal and also for the continuum and 
b ^ c backgrounds. The fit result is ^obs = 0.05 ± 0.16, 
corresponding to a w-corrected asymmetry Acp = 0.06± 
0.19. 

The charge symmetry of the detector performance and 
reconstruction procedure is verified with a sample of 
_B+ -^ £'°7r+, Z?" -^ K^Ti~'iT^ decays and their charge 
conjugates. We apply the same procedure that is used 
for 5+ — > p+7r° to select B^ -^ D^n'^ candidates and 
extract signal yields by fitting the /S.E distribution. For 
these events we find a direct-CP-violating asymmetry of 
—0.03 ± 0.02. We assign 0.03 as the systematic error as- 
sociated with detector and reconstruction effects. The 
systematic error associated with the fitting procedure is 
determined to be 0.01 by shifting each fitting function 
parameter by ±lcr and taking the quadratic sum of the 
resulting changes in Acp- An error of lo;o5 fo'^ t^oti- 
resonant background is estimated by subtracting the non- 
resonant component, obtained from the fit to M(7r~7r°) 
or M(7r+7r°) and increased by ±lcr, from the B~ or B^ 
signal yields. The change in Acp is taken as the system- 
atic error. A 0.01 systematic error due to the background 
from other rare B decays is assigned. The quadratic sum 
of these errors is taken as the total systematic error. We 
obtain a CP-violating charge asymmetry that is consis- 
tent with zero: 



AcpiB-^ 



p^n") 



= 0.06±0.19(stat.)to:°^(sys. 



In summary, we observe the decay B^ -^ p+vr" with a 
statistical significance of S.lcr and measure its branching 
fraction. The results are consistent with those from 
the BaBar experiment |9(. This measurement provides 
one of the essential quantities to constrain 02 from an 
isospin analysis oi B ^ pn decays. The measured charge 



asymmetry is consistent with zero. 
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